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Abstract 

We calculate a two-loop effective potential to the order of 0{X^as) in the 
MSSM. We then study the corresponding two-loop corrections to the CP-even 
Higgs-boson mass for arbitrary tan/3 and left-right top-squark mixings. We 
find that the lightest Higgs-boson mass is changed by at most a few GeV. We 
also show the improved scale dependence and compare to previous two-loop 
analyses. 
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1 Introduction 


In the minimal snpersymmetric standard model (MSSM), the Higgs sector is com¬ 
posed of three nentral (two CP-even, one CP-odd) and two charged scalar bosons. 
An important fact of the model is that the qnartic self-conpling of the lightest CP- 
even Higgs boson h is not a free parameter, bnt related to the standard model gauge 
couplings g and g'] as a result, the Higgs boson has an upper bounded tree-level mass, 

< Mz- This limit, however, is violated when the one-loop radiative corrections 
are included. 

The well-known dominant one-loop radiative correction comes from the incom¬ 
plete cancellation of the virtual top-quark and top-squark loops |^, it approximately 
has the size Am| ~ log(m?/m^). Taking = 100 — 1000 GeV, one hnds 

a large correction. Amt — a few — 50 GeV, due to the relatively large top-quark 
mass = 175 GeV. The one-loop Higgs-boson mass sensitively depends on the 
top-quark mass, and generally varies with the change of renormalization scale. So 
it remains a quite important problem to study the magnitude of two-loop radiative 
corrections and the scale dependence of the Higgs-boson mass after these corrections. 

There are basically two approaches to calculate the two-loop radiative corrections. 
In the renormalization group equation (RGE)-improved effective potential approach 
iii. one uses the one-loop effective potential, together with two-loop RGBs, then 
all leading-order and next-to-leading-order corrections can be calculated. The hnite 
one-loop threshold corrections arising from the decouplings of the heavy top-squarks 
have also been included, both for the small and large left-right top-squark mixings. 
It is further observed that by judiciously setting the renormalization scale, one can 
use the one-loop renormalized Higgs-boson mass as a good approximation to the full 
two-loop results ia- 

The second approach involves a two-loop effective potential. In the special case of 
tan /3 —> oo and no left-right mixing in the top-squark sector. Hempfling and Hoang 
have calculated the upper bound to the two-loop Higgs-boson mass [|] . Their results 
qualitatively agree with the previous approach. Recently, two-loop corrections to ruh 
have also been computed by an explicit diagrammatic method 0]. 

It is the purpose of this paper to generalize the two-loop effective potential cal¬ 
culation of [^] to the case of arbitrary tan/? and left-right top-squark mixings. The 
effective potential method has an advantage over other approaches because it is sim¬ 
ple and does not require complicate programming. In this paper we will study the 
improved scale dependence of the Higgs-boson mass rrih and the size of two-loop 
corrections. We will also compare our results with previous two-loop calculations 
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The rest of the paper is organized as follows: We hrst present a general formalism 
for calculating the GP-even Higgs-boson mass from the effective potential in Section 
2, and compute the two-loop effective potential to the order of O^Xfus). We next 
show the results of our numerical analyses in Section 3; we hnd good agreements with 
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previous two-loop calculations. Finally we conclude in Section 4. For completeness, 
some functions which appear in the two-loop calculation are given in the Appendix. 


2 Effective Potential and the CP-even Higgs- 
boson Masses 


We start our analysis with the tree-level potential of the MSSM^, 


Vtree — + fX^)\H2\‘^ + fiB {H1H2 + 

,2 


+ g'^ 


(Pi 


\H, 


| 2\2 




( 1 ) 


where g,g' are the SU{2) and f/(l)y gauge couplings, mHi,rnH 2 and B are the soft- 
breaking Higgs-sector mass parameters, and p is the supersymmetric Higgs-boson 
mass parameter (the /i-parameter). 

We express Hi and H 2 in terms of their component fields. 


Hi = 


{Si+iPi)/^ 

Hi 


Ho = 


Ht 

{S2 + lP2)/^ 


( 2 ) 


so the tree-level potential can be rewritten as a function of the CP-even fields. Si 
and 82 - In general, the all-loop effective potential is a function of Si and 82 , which 
are usually known as classical fields. 

The technique for calculating a higher loop effective potential was developed long 
ago by Jackiw 0. First, the Higgs fields are expanded around the classical fields, in 
terms of which all the relevant particle masses and couplings are determined. One 
then calculates the higher loop effective potential by computing the corresponding 
zero-point function Feynman diagrams (bubble diagrams). 

To be more specific, to the two-loop order that we consider in this paper, we write 
the effective potential as 


1 ^( 51 , S'2) — Vo -|- ldiree(*S'l, 52) -|- Fiioop(*S'i, S'2) -|- V 21 oop(>S'l, 52) , 


( 3 ) 


where Vq is a field-independent vacuum-energy term, and Vtree, Viioop and V 2100 P are 
the tree-level, one- and two-loop contributions respectively. Vq is necessary for the 
renormalization group invariance of the effective potential . 

The one-loop effective potential in the DR scheme is well-known. It can be easily 
obtained by calculating the one-loop bubble diagrams with all kinds of (s)particles 
running in the loop. We find that in the Landau gauge 


ri,.„p(Si,S2) = Y.Y.NlGCfi}-2Y.NlG(f) + 3G(W) + -G(Z) 

/ i.i / ^ 


^We work in the modified DR -scheme of Ref. 
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G{H) + G{h) + G{A) + G{G) 


+ G{H+) + G{G^) 


2j:G{xt)-T.Gm > 


( 4 ) 


2=1 


2=1 


where / sums over all the (s)quarks and (s)leptons, iV/ is the color factor, 3 for 
(s)quarks and 1 for (s)leptons, and all the masses are implicitly Si, S' 2 -dependent. H, 
h, A and G and G'*') label the neutral (charged) Higgs and Goldstone bosons, 
Xi and represent charginos and neutralinos. We have also used short-handed 
notations fi = , W = etc. The function G{x) is dehned as 

Ji 

where In a: = ln(x/(5^), with Q the renormalization scale. 



Figure 1: Bubble diagrams for the two-loop effective potential to the order of 0{X^as) 
in the MSSM. 


The two-loop effective potential can be derived similarly, the corresponding Feyn¬ 
man diagrams are plot in Fig. |^. To the order of Oi^X'fas), we have in the Landau 
gauge 

V2ioop(5'i, S 2 ) = 327ra^| - 2t l{t,t,0) 

1 

+ -(c^ -f- s^) ti) -I- 2sjC^J(fi, ^ 2 ) + ti, 0) 

^ i=l i=l 

+ J2L{ii,g,t) -Am-g rrit StCt(^I{ii,g,t) -/(f2,^,t))| (6) 

where the (minimally) subtracted functions I and J are dehned in the Appendix, st is 
the top-squark mixing angle. This effective potential, in the limit of tan > cx) and 
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no left-right squark-mixing (st = 0), agrees with that of Ref. 0. As a good check, one 
can show that the effective potential ^(Ri, 5'2) is invariant under the renormalization 
scale change, up to two-loop terms which are ignored in our approximation. 

The effective potential Eq. (P), as a generating functional, encodes the information 
of two-loop tadpoles and self-energies at zero external momentum. From this, one 
can End two-loop CP-even Higgs-boson masses by solving appropriate on-mass-shell 
conditions. 

This proceeds as follows: We hrst minimize 1/(51, 52 ) at the Higgs vacuum ex¬ 
pectation values Si = Vi and S 2 = V 2 , 


dV 

Wi 


Si=vi,82=^2 


0, 


dV 

W 2 


Si =V1,S2=V2 


0 , 


which is equivalent to the requirement that the tadpoles vanish. 


(7) 


Ti 1 

— = -m%C 2 i 3 + + 1 /+ Bfi tan f3 , (8) 

T 1 

— = -o"^|c2/3 + -h P/i cot 13 , (9) 

Vo Z 


where tan/3 = V 2 /V 1 , m| = + g'‘^)v‘^/4: and + v\. To the two-loop order, 

the tadpoles T*, z = 1, 2 are given by Ti = 7 t 1 '°°p_|_ 2721 oop^ where the one- and two-loop 
tadpoles are dehned by 


T 


Hoop 


91/lloop \ 

dSi J 


Si=i;i,S2=^'2 


T 


21oop 


^lAloop \ 

dS, J 


Sl=Vi,S2=V2 


( 10 ) 


One can check that the one-loop tadpoles obtained in this way give the same 

results as in Ref. OB where they were explicitly calculated from the one-point 
function Feynman diagrams (tadpole diagrams). 

The CP-even Higgs-boson mass matrix, after some algebra, is 

M^(p^) = 

m%cl - Renii(p2) + Ti/vi -(m| m\)sfsCp - Reni2(p^) 

-(m| m\)sf}C 0 - Reni 2 (p^) m|s| -F - Ren 22 (p^) + T 2 /V 2 

where m\ = —pP(tan (3 -|- cot /3), and both mz and are DR running masses. H’s 
represent two-point functions (self-energies). 

The radiatively corrected Higgs-boson masses can be found by computing the 
zeroes of the inverse propagator, The complete one-loop formulae for 

Hn Ref. |]ll| , we used the ’t Hooft-Feynman gauge. To compare with the one-loop tadpoles 
obtained from Eq. (^), the Goldstone boson masses need to be set to zero. 
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self-energies at nonzero external momentum can be found, e.g, in Ref. 
two-loop self-energies, we shall approximate them as follows: 


111. For the 



fd^V21oop\ 

V dSidSj J 


1 

<S'i=Ri,52=R2 


hi = 1,2. 


( 12 ) 


The difference of and (p^) is negligible. 


3 Numerical Procedure and Results 


We shall hrst show the improvement of renormalization scale (Q) dependence of the 
lightest CP-even Higgs-boson mass. We start by solving the two-loop renormalization 
group equations (RGBs) of the MSSM. The boundary condition inputs for these RGBs 
are taken to be the observables aem, Gp, Mz, mt, mb and m,- at the electro- 
weak scale, Mz, and the universal scalar soft mass Mq, gaugino soft mass M 1/2 and 
trilinear scalar coupling Aq at the unihcation scale, Mqpt ^ 2 x 10 ^® GeV (This 
is sometimes called the minimal supergravity (mSUGRA) scenario.). We have also 
properly taken into account the low energy threshold corrections to the gauge and 
Yukawa couplings, as in Ref. Q- 

In this scenario, the //-parameter and the (tree-level) mass are determined 
in terms of other variables through the minimization conditions (|^) and (^, which 
ensure the correct electro-weak symmetry breaking. 


= 


m\ = 


tan 2/3 ( P ~ cot {3 


Hi 


— mr 


cos 2/5 


—2 —2 \ 2 

-^Hi) -^Z , 


(13) 

(14) 


where = mjj^ —Ti/vi and = m\j^ —T 2 /v 2 - 

We then calculate the lightest CP-even Higgs-boson mass from the mass-squared 
matrix Bq. (0). We use the one-loop tadpole and self-energy formulae from Ref. 0], 
while for the two-loop contributions, we compute the tadpoles and self-energies 
numerically according to Bqs. (0 and ( 0 ) by replacing the differentiation by a 
hnite difference. The held-dependent masses in Bq. (^ are the top-quark mass 
mt = XtS 2/\/2 and the top-squark masses which are found from the following held- 
dependent mass-squared matrix: 


/ + 1A?S| + ^Si) \ 

( 75A,(3l,S2 + mSi) Mfj + j 


where we have neglected the held-dependent P-term contributions and Mq, Mp are 
squark soft masses at the low-energy scale. The held-dependent angle St in Bq. (P) 
is dehned as the mixing angle of the above mass-squared matrix. 
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Figure 2: Renormalization-scale (Q) dependence of the lightest GP-even Higgs-boson 
mass rrih. The dashed and solid lines correspond to the one- and two-loop masses 
respectively. We have hxed the universal boundary conditions, Mq = 500 GeV, 
Mi /2 = 200 GeV, Aq = 0, and chosen a negative p-parameter. 


In Fig. 1^, we show the dependence of one- and two-loop radiatively corrected GP- 
even Higgs-boson masses mh on the renormalization scale Q. We choose the universal 
soft parameters Mq = 500 GeV, M 1/2 = 200 GeV and Hq = 0 at the unihcation scale, 
and set the sign of p-parameter to be negative^]. We plot two choices of tan/3, 2 and 
20. The one- and two-loop masses are shown in dashed and solid lines respectively. 
In the corrections to the Higgs-boson masses nih, we have used the DR running top- 
Yukawa coupling at the scale Q. The formulae which convert the top-quark pole 
mass to Xt{Q) are given in Ref. [0- 

We see that the one-loop radiatively corrected Higgs-boson masses vary by about 
10 GeV as the renormalization scale Q varies from 100 to 500 GeV. However, once 
we properly include the two-loop radiative corrections, the scale dependence of rrih 
becomes much milder, and we hnd for all ranges of the scale rrih — 88 GeV and 110 
GeV for tan /3 = 2 and 20 respectively. The two-loop calculation can only change mh 
by a few GeV. 

In Figs. ^ and |^, we compare our results with that of RGE-improved effective 
potential (EP) approach IQ. In that approach, the heavy particles are decoupled at 
m/ (or m/^ and stepwisely if they are very different). The two-loop RGEs of 
the effective held theory below the decoupling scale is then used to run the effective 
couplings to the scale where the Higgs-boson mass is evaluated, e.g. the on-shell top- 

^Our convention of the sign of ^-parameter is opposite to that of Refs. 0 and 0. 
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Figure 3: Higgs boson masses rrih vs. the squark soft masses Mg in the no-squark- 
mixing case = /i = 0. The solid lines are results from the two-loop EP approach. 
For comparison, we also show the results from the RGE-improved one-loop EP ap¬ 
proach in dashed lines. 


quark mass. Since the next-to-leading-order corrections are negligible at this scale 
PI, this allows an analytical solution to the two-loop RGEs 0. 

In this part of numerical analysis, we do not impose the minimization conditions 
Eqs. (ID and (|D, instead, the /i-parameter and the GP-odd Higgs boson mass tra 
are taken as inputs. We further choose the squark soft masses Mq = Mu = Mg. 
In Fig. 1^, we show the two-loop Higgs-boson masses mh versus Mg, in the RGE- 
improved one-loop EP approach (dashed lines) and in the two-loop EP approach 
(solid lines). The parameters and p are set to zero; this corresponds to the no 
left-right squark-mixing limit. Other parameters are = M 3 = Mg. Results 
from the two approaches generally agree to < 3 GeV, for both small and large tan (3 
cases. The difference comes from the 0{Xf) corrections which are neglected in this 
calculation p. 

In Fig. we choose a nonzero /i-parameter, /i = —200 GeV, and plot rrih versus 
Xt/Mg, where Xt = At + tan/? is related to the off-diagonal element of the squark- 
mass matrix. We see the results from the two approaches still agree quite well except 
for the large Xt/Mg. In particular, the curves for the two-loop EP approach peak 
at Xt/Mg ~ ±2, which is different from the RGE-improved one-loop EP approach. 
This however agrees with the results of a recent analysis 0]. The slight asymmetry 
of the solid curves with respect to the reflection, Xt —> —Xt, originates from the last 
term in Eq. (|D- We found that for tan/? = 40 the upper limit for the Higgs-boson 
mass is 125 GeV at the large squark-mixing [Xt/Mg = ±2). 





Xt/Ms 


Figure 4: Higgs-boson masses rrih vs. Xt/Ms, where Xt = At + n/ tan/3. Both the 
results from the RGE-improved one-loop EP and two-loop EP approaches are shown. 


4 Conclusions 

To conclude, we have used an effective potential method to calculate the two-loop 
corrections to the lightest GP-even Higgs-boson mass in the MSSM. Our approach 
is straightforward and easy to program, and can be extended to include two-loop 
corrections of order 0{X/). We show that the renormalization scale dependence 
of rrih improves after including the two-loop corrections, this largely reduces the 
uncertainty associated with one-loop calculations. We have shown that the two-loop 
correction is only about a few GeV with respect to the one-loop results (where the 
DR running coupling Xt is used). We have also compared our results with some 
previous two-loop calculations. We found good agreements with the RGE-improved 
one-loop EP approach of Ref. @] except for the case of large left-right squark mixing, 
where we obtained similar results as in Ref. 0 . The upper bound for the Higgs-boson 
mass nih ^ 125 GeV is achieved at the region of parameter space for large tan f3 and 
left-right squark mixings. 
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Appendix: The functions I{x,y,z) and J{x^y) 


Momentum integrals arising from the two-loop bubble diagrams have one-loop 
sub divergences which can be subtracted in the standard way Furthermore, in 
the DR scheme there is no complication associated with vector-boson loop^ , so all 
integrals can be expressed in terms of (minimally) subtracted functions J, J which 
are 


(levr ) J{x,y) = X y 


1 — In a: — In 1 / -|- Inx In y 


(A.l) 


(167r ) I{x,y,z) = -- 


{y + z — a;)ln y In z + {z -\r x — y)\n. z In x (A.2) 

^ L 

-|-(x + y — z)\\i X hi y — A{x In X + y \n. y + z In z) + ^{x, y, z) + 5(a; + y + z) , 
where ^ is given by 

^{x,y,z) = 8b 

when —6^ = a? = {x"^ + y"^ + z"^ — 2xy — 2xz — 2yz)/4 < 0, and 


71 


L{e,) + LiOy) + L{e,) --\n2 


(A.3) 


^{x,y,z) = 8a 




when > 0. Here L{t) is Lobachevsky’s function, dehned as 

1^, ,;i,_iSin2H 


Pt I ^ 

L{t) = — dx In cos a; = tln2 —-^( 

2 k=i 

and the function M{t) is dehned as 




M(t) = - f 

an 


dx In sinh x = 


TT 

12 


^+ tln2 - ^Li 2 (e 


(A.4) 


(A.5) 


(A.6) 


Li 2 is the dilogarithm function. 

The angles and (l)x,y,z are dehned by 

'y + z — X 


9r = arctan 


26 


= arccoth 


y + z — X 
2 a 


etc. 


(A.7) 


Finally we have also used the following function in the two-loop ehective potential: 
L{x,y,z) = J{y,z) - J{x,y) - J{x,z) - {x - y - z)I{x,y,z). (A.8) 


^In contrast, the vector bosons live in d = 4 — 2e dimensions in the MS scheme, in reducing the 
two-loop integrals to the form of / and J, there will be extra finite terms from themselves and the 
associated one-loop subdiagrams. 
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